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SUMMARY 

Matrix methods are employed in the solution of absorption 
problems In vhich the Incident radiation is a known arbitrary 
function of time. The case of a plane source of polychromatic 
radiation of several types at normal incidence to a plane absorber 
is considered. 

The method Is applied to a hypothetical example in vhlch the 
amounts of radiations transmitted and reflected as veil as the 
amounts converted to thermal energy are found for a time -dependent 
source. 


INTRODUCTION 

A mathematical analysis of the interaction of nuclear 
radiations vlth matter has been applied to the case where radio* 
activity of the absorber is time-dependent (reference 1) . The 
intensities of sources of nuclear radiations, however, may also 

be time -dependent. For example, U 235 , an alpha emitter vlth a 

8 231 

half-life of 7.1 X 10 years decays exponentially to Th , 

a beta emitter with a half-life of 24.6 hours that in turn becomes 

Pa 231 , another alpha emitter vith a half-life of 3.2 x 10 4 years. 
By the time the decay process is complete, the final product 

being stable Pb 207 , the original U 235 has become a mixture of 
13 isotopes of 11 elements and the nuclear radiations that are 
emitted from this mixture consist of alpha rays, beta rays, and 
gamma rays of various energies and intensities. Thus, an absorber 
designed to protect personnel from the effects of the nuclear 
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radiations emitted by U 23 would be inadequate unless it could 
adequately diminish all these radiations and also the radiations 
due to radioactivity induced in the absorber by all the radiations. 

58 

Conversely, Cu , a positron emitter with a half-life of 

81 seconds, decays exponentially to stable Ni 58 and consequently 
would require no shielding whatever after 1 day - an obvious 
result that is not predicted in a theory that assumes constancy of 
source intensity. 

In order to enable shielding calculations of such time- 
dependent sources to be made, an extension of the methods of 
reference 1 to include the case in which the incident radiation 
is a known arbitrary function of time was developed at the NACA 
Lewis laboratory and is presented. Explicit use of the more 
general case, in which the radiation absorbed may -consist of any 
finite number of types, is shown herein. 


SYMBOLS 


For the most part, the definitions of all the symbols used 
agree with the definitions in reference 1. Brackets around the 
matrix terms have, however, been omitted whenever no decrease in 
clarity results. The following symbols are used in this report: 

[A],[B],[c],[D] matrix coefficients 

F(t ) matrix whose elements are known functions of time 


H 

h 


li 

✓ 


n 

P 


R 


i 


matrix of thermal power generated 

matrix of rate of conversion of energy of nuclear 
radiation absorbed to thermal energy 

matrix of power of radiation from (i-l) th station 
incident on i - ^ station 

number of stations in absorber 

matrix whose elements are one -half of power of 
radiation from radioactivity 

matrix of power of radiation from i - ^ station 
incident on (i-l)^ station 

matrix of power back- scattering coefficient 


r 
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S i 

t 

X 

»*i 

T 

Subcripts : 

0 

i 

T 

Superscripts : 

fi • f • • • $ t • • 

n 

7 


i^ 1 station of at sorter 

matrix of power-transmission coefficient 

matrix of decay constant for radioactivity, 
second • 

matrix of rate of conversion of energy of 
radiation at sorted in Sj to energy of 
radioactivity in S^ 

time, seconds 

initial condition 

i^* 1 station of absorber 

time 

types of radiation 
neutron radiation 
gamma radiation 


ANALYSIS 

The analysis is tased on the following assumptions: 

(1) Elements of absorber are parallel plane surfaces. 

(2) Incident radiation is a known arbitrary function 
of time . 

(3) Incident radiation consists of several types of poly- 
chromatic radiation. 

(4) Radiation is normal to atsorter. 

(5) Part of radiation absorbed at each station of absorber 
is transformed to thermal energy and part to induced radioactivity 
in the absorber. 
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(6) Badioactive atoms in the absorber decay in a single step 
to stable atoms. 

(7) Products of radioactive decay within the absorber may 
consist of types and energies that are the same or different from 
those of the incident radiation. 

(8) One-half of the radioactivity produced in each station 
of absorber is emitted from each side of the absorber. 

(9) Radiation impinging upon source R-j_ has a negligible 
effect on the source intensity. 

Method . - Diagram (a) is included to clarify the discussion 
in the following paragraphs. 

S i 

The absorber may be considered as 
consisting of n stations normal to the 
path of radiation. The i^h station of the 
absorber is denoted by Sj_. Radiation I± 
from the preceding station Sj_ _i is 
incident upon one surface of whereas 

radiation Rj+i from the subsequent 
element S^ + i impinges upon the opposite 
surface of S^. Prom Si emerges radiation 
Ii+1> ' wh ick is incident upon and 

radiation Ri, which is incident upon H i 

(a) 


In the treatment presented herein, the following quantities are 

assumed to be known: 

I-]_(t) radiation power incident upon absorber from source 

]\i+i(t) radiation power incident upon absorber from side of 

absorber opposite source 

(Pi)o one-half of initial power of radioactivity in i^ h station 
of absorber 

(H-^q initial thermal power being emitted from i^h station of 
absorber 
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h matrix of rate of conversion of energy of nuclear 

radiation absorbed to thermal energy 

r matrix of power back-scattering coefficient 

t matrix of power -transmission coefficient 

X matrix of radioactive decay constant 

U matrix of rate of conversion of energy of radiation 

absorbed to energy of radioactivity 

All of the previous quantities except the first two are 
time -independent constants. 


The problem is to find 1 ^ , and in terms of known 

quantities. The defining equations and conditions for q types 
of radiation are: 


ii + i a ■ *i a ii a + r i a ®i + i a ♦ V 

B i “ r i x i + *1 B i+i + p i 
dPi 1 * * 




_i_ = (I* + E 1+1 a ) ♦ . . . ♦ <V + 8 1+1 J ) + • 




h+l* " t i J V + r i J B i + l J + V 

B i J “ r i J V + t i J E i + l J + P i J 




(1) 


(Ij 4 + E 1+1 a ) + . . . + + R 1+1 J ) + • 

dl 1 + Bi + l 4 ) - X i J Pi* 5 
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*U1 “ *1* x i 4 + r i q B i + i* * p i q 


= r * I * + t * R 


i+1 


+ P, 


^ , *'* (I ± * + R 1+1 a ) ♦ 


+ H. 




(v + R i + 1 } + • 


V (l) Continued 


Hi*' 4 (I^ + B 1+1 4) - Xj4 Pi 4 

Hj ■ h^ a + ®i+i*) + . . . + (Ij^ + ®i+i^) + • • . + hj 4 (I^ 4 + Rj + i 4 ) 

J 

where the double superscripts signify a conversion of radiation of 
the type represented by the first superscript to the type repre- 
sented by the second superscript. 


i 

Equations (l) are used to find P. As in Case HI of refer- 
ence 1, all the I's and E' s are expressed in terms of I^(t) 
and I^i+i(t), which are known functions of time and of P's. 

The results may be expressed in the form 


. w &] . [b] [ijcol + [c] i^ + 1 (t) + ta (z) 


where [ A] , [ij] f [c] , and [d] are matrices, all the elements of which 
are constants. 


Let 

[Bljl-L^jj + [C][% +1 (T)J + [D] =^F(t)J 


Then if [Ii(t)] and [^^(t)] are independent of [p] , the 
solution to equation (2) may be written as 



(3) 


Thus, [p] is reduced to a quadrature. Substitution of the 
expressions for [jl( T )] ©*“* |®n+l( T J ©Ration (3) results 

in an explicit expression for P. After P is known, Ij_, Ej, and 


Ej_ may be found. 
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Previous knowledge of all the nuclear reactions that will 
occur during the time interval in which the calculations are made 
is necessary; for if even one energy or type of radiation is 
present hut not accounted for, all the calculations may he 
invalidated. 


Examgle. - Assume that: 

(1) The absorber is two stations in thickness. 

(2) Energy is not degraded nor upgraded. 

(3) The incident radiation .is 1.00 Mev gamma radiation. 

(4) The Incident gamma radiation induces a radioactivity of 
2 Mev neutrons in the second station of the ahsorher. 


(5) I x 


0) 



hi y»y 




pirl 

' s 




^ ® in WT ) roentgen per hour 


ttx 10~ 3 second -1 

roentgen per hour 


= 0.5 


= 0.8 
s 0.9 
= 0.9 
= 0.1 
= 0.05 


- 0.1 


= 0.1 

= \i-2jf >7 = Hi?* 11 = (ij 11 # 11 = ~ M2 n * n 


0 
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[x 2 7,n = 10"® second” 

X 2 n = 10"® second” 1 

h l 7 ■ h i n ■ h 2 n ■ 0 
• h 2 7 * 10" 3 

Mo* Ho- Ho-° 

At T= 0, P 2 n = 0.1 roentgen per hour. The problem is to deteimine 
the values of I 2 , T 3 , R 1# R 2 , and H 2 at T = 10 5 seconds. 

Solution . - If equations (l) are applied to the previous data 
and brackets are emitted for simplicity, then 

I 2 7 > 0.5 I-J + 0.1 H 2 7 

I 3 7 =0.8 I 2 7 

= 0.1 I-J + 0.5 R 2 ? 

R 2 7 =0.05 I2 7 

or in terms of I] 7 , 

I 2 7 = 0.503 Ij 7 
I 3 7 = 0.402 I x 7 
R x 7 = 0.113 I-, 7 
R 2 7 = 0.025 Ij 7 

Also 

I 2 n =0.1R 2 n 
I 3 n = 0.9 I 2 n + P 2 n 
R^ * 0.9 R 2 n 
R 2 n = 0.1 I 2 n + P 2 n 



from which 


I 2 n = 0.101 P 2 n 

I 3 “ = 1.091 P 2 n 

= 0.909 P 2 ? 
R 2 n = 1.010 P 2 n 


Now 


dP 2 n ft 1 ^ c . 

- — = 10 b second " 1 x 0.503 Ij? - 10 " 5 second " 1 P 2 n 


Or 


P 2 “ = e- 10 - 5T < 


0.1 + Is, 05 x 10 " 7 e “ 70 Sr I 1 l(T) 


dT 


= e 


■l°" 5r <|o.l + 5.03 x 10" 7 /C 10 “ 5t (1 - ain cot) dT 


q -10-5t _ 

— Jq + 5.03 x 10"' 


1 0) cos (i) t - 10~ 5 sin 

10- 5 c 2 . 10 -10 


At t =10 seconds. 


P 2 = 0.0873 roentgen per hour 


Therefore 


- 0.503 roentgen per hour 

y 

1*3 = 0*402 roentgen per hour 

y 

R^ = 0.113 roentgen per hour 
Rg^ = 0.025 roentgen per hour 
I 2 n = 0.0088 roentgen per hour 
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0.0952 roentgen per hour 
0.0794 roentgen per hour 
0.0882 roentgen per hour 
10“^ = 5.03 x 10 roentgen per hour 


DISCUSSION 

In practice , most sources of radioactivity are time -dependent. 
This time -dependency varies from the simple exponential natural 
decay law obeyed by all radioactive substances to the sharply vary- 
ing functions of dynamic piles. Thus, there are many cases in 
which the amount of shielding material may safely be less than 
that required to diminish adequately the raavi 7m m instantaneous 
output of a source. Application of the method described herein 
provides a frequently quick answer to the amount of shielding 
needed for a given source. 

In general, the right side of equation (3) includes an 
integral of an arbitrary function. In such cases, approximate 
and numerical methods of integration may be applied to obtain 
explicit values for P^. This method requires no further 

integration and all values of I}, and % can be found by 
elementary algebraic manipulations involving . 

The main limitations of the method are a lack of sufficient 
experimental data on the absorption of products of radioactivity 
and on the rates of conversion of these radiations to secondary 
radiations and to thermal energy as in reference (l). Consider- 
ation of the geometry must also be made; a plane source at normal 
incidence to a plane absorber is assumed throughout. The magnitude 
of the effects of deviations from this ideal case must be known 
for reliable results in a given set of calculations. 


SUMMARY OF RESUITS 

Matrix methods are employed in the solution of absorber 
problems in which the incident radiation is an arbitrary known 
function of time. The solution is obtained in the form of a 
quadrature. 



MCA TN 1952 


11 


The method presented is applied to a hypothetical example in 
which the incident radiation consists of a source of g^rrmw. radiation 
that varies sinusoidally with the time. 
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